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Abstract 

g 

Several realistic phenomenological nucleon-nucleon interaction models are employed 
[ — ■ to investigate the possibility of bound deuteron-like states of such heavy flavor 

hyperons and nucleons, for which the interaction between the light flavor quark 
I components is expected to be the most significant interaction. The results indicate 

that deuteron-like bound states are likely to form between nucleons and the H c and 
H cc charm hyperons as well as between E hyperons and double-charm hyperons. 
Bound states between two S c hyperons are also likely. In the case of beauty hyperons 



the corresponding states are likely to be deeply bound. 



o 



1 Introduction 



O 

The interaction between nucleons and several classes of heavy flavor hyperons 
is expected to be dominated by the long range interaction between the light 
flavor components of the baryons. As an example, the interaction between 
nucleons and the recently discovered double-charm hyperons [1,2,3] is mainly 
due to the interaction between the single light flavor quark in the double fla- 
vor hyperon and those in the nucleon. The color-neutral interaction between 
charm and light flavor quarks is either weak or of very short range. In the case 
of two-baryon states, in which the interaction between the light flavor quarks 
is the dominant one, it should be possible to calculate the binding energy to a 
first approximation by modifying realistic phenomenological nucleon-nucleon 
interaction models to take into account the different numbers of light flavor 
quarks in the baryons [4] . This method is explored here. It is found to have at 
most qualitative value, however, due to the fact that the short range compo- 
nents of the extant realistic phenomenological nucleon-nucleon potentials are 
poorly constrained. Its use should be limited to those potentials which do not 
have strong angular momentum dependence. 
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The modification of the nucleon-nucleon interaction required for systems with 
different numbers of light flavor quarks is straightforward when the nucleon- 
nucleon interaction is expressed in terms of operators, which have well defined 
matrix elements in the quark model, and the strength of which may there- 
fore be correspondingly rescaled from the nucleon-nucleon system to the two- 
baryon system under consideration. In this approach the radial behavior of 
the interaction components is determined by the phenomenological nucleon- 
nucleon interaction, without any need for a detailed microscopic quark model 
based derivation of the interaction. Two-baryon systems that may — at least 
approximately — be described in this way are the following: iV — E c , N — S c , 
iV — S cc as well as those with the nucleon (N) replaced by the paired charm 
hyperon and finally also the corresponding states that involve beauty rather 
than charm hyperons. To the extent that the interaction between the strange 
and light flavor quarks may be neglected the interaction between S hyperons 
and the corresponding charm hyperons may also be described in this way in 
a first approximation. In contrast the interactions of heavy flavor hyperons 
with zero isospin like the A c and the A 6 hyperons, which do not couple to 
pions, cannot be approximated by rescaled versions of the nucleon-nucleon 
interaction. 



The interaction between nucleons and charm and double-charm hyperons is 
weaker than that between the nucleons, because of the smaller number of light 
flavor quarks in the charm and double-charm hyperons. The weaker attractive 
interaction is largely compensated, however, by the weakening of the repulsive 
effect of the two-baryon kinetic energy that is caused by the larger masses 
of the charm and double-charm hyperons. As a consequence the calculated 
binding energies remain small and dependent on the details of the interaction 
model in the case of charm hyperons, and typically become large only in the 
case of systems of beauty hyperons. 



Here the interaction between the two-baryon states that are formed of nucleons 
and heavy flavor hyperons is calculated from the realistic phenomenological 
interaction models in Refs. [5,6]. For an estimate of the theoretical uncertainty 
in the calculated binding energies the set "AVn" of systematically simplified 
versions of the AV18 interaction model is also employed [7] . 



Section 2 below contains a description two-baryon states that are formed of 
nucleons and iso-doublet hyperons. Two-baryon states with isospin 1 hyperons 
are considered in section 3. The results are summarized in the concluding 
discussion. 
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2 Two-baryon states of isospin 1/2 baryons 

Two-particle states formed of heavy flavor isodoublet baryons are similar to the 
two-nucleon system in that they have a long range pion exchange interaction, 
and an attractive intermediate range component, which to a large extent may 
be attributed to two-pion exchange [8,9]. To the extent that the color-neutral 
interactions between their heavy flavor quarks is weak in comparison to the 
interaction between the light flavor quark components, the interaction may - 
to a first approximation — be constructed by multiplication of the components 
of the nucleon-nucleon interaction by appropriate quark model scaling factors. 

The phenomenological nucleon-nucleon interaction is in general expressible in 
terms of rotational invariants of spin and isospin operators as well as momenta 
and angular momenta. The scaling factors for the strengths of the matrix 
elements of spin and isospin invariants for two baryon states of isodoublet 
baryons relative to the corresponding two-nucleon state matrix elements may 
be derived from the quark model matrix elements of the spin and isospin 
operators for light flavor quarks. For S cc (and S) hyperon states the matrix 
elements are: 



( S (cc) 1 


E 1 ) = 

q 


\(N\^l q \N), 

6 q 


( S (cc) | 


E°f 1 ) = 

q 


~(N\^\N), 
6 q 


( S (cc) 1 


E T i 1 E (cc) ) = 

q 


q 


( S (cc) 1 


E a " r l 1 s (-) ) 

q 


= ~(N\^7j\N) 
q 



For the corresponding single charm hyperon states formed of Sj, hyperons the 
matrix elements are: 

(Z' c \^ q \Z' c ) = \(N\Y,l q \N) , 

q q 

(s;ie^is;> = |we^iao, 

q <3 q 

(K\Y,t?\K) = (n\Y,t?\n), 

q q 
(K\E^\K)^1(N\J:^tJ\N). (2) 

q q 

From these matrix elements one may derive the scaling relations that apply 
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for the rotational operators in nucleon-nucleon interactions. As an example 
the AV18 interaction model is given in the operator form [6]: 

18 

V = J2Mr)OV; (3) 
p=i 

OP =i,..i4 = [i^. . ^, ^, L • S, L 2 , Pa, ■ d v (L ■ S) 2 ] ® [1, n ■ r 3 ] . (4) 

The operators O p=15 '~ 18 are part of the isospin breaking electromagnetic com- 
ponents which depend on the baryon charge. These (small) terms are dropped 
here since they are inapplicable to hyperons. 

The decomposition used in Ref. [5] is slightly different from (4): 

p=i,.w = [x^. . ^, Sij, L ■ S, Qij] ® [1, n ■ tj\ , (5) 

where = \[{(Ji ■ L){Bj ■ L) + (a,- • L)(ai ■ L)]. The corresponding scaling 
factors are listed in Table I. 

The total spin operator S = ((7j + <7,-)/2 does not yield an unambiguous scaling 
factor, when the matrix elements contain terms with different scaling behavior. 
In the first and the third column of Table 1 (bound states with nucleons), 

— * 

the scaling factors for all operators, which contain the total spin-operator S, 
have been approximated by the largest scaling factor in the corresponding 
expression. In the other columns (bound states of two hyperons) only the 
scaling factors containing (L ■ S) 2 had to be approximated in the same way. 

The scaling factors in Table 1 reveal that the interaction between E,' c and S cc 
hyperons differs qualitatively from the nucleon-nucleon interaction in that the 
spin-independent central interaction, which contains most of the intermediate 
range attraction, and the strong short range repulsion is weaker by an order 
of magnitude. This weakening of the short range repulsion increases the rel- 
ative importance of the other short range interaction components, which are 
of little significance for — and as a consequence not very well constrained 
by — low energy nucleon-nucleon scattering observables. As these other short 
range components vary significantly between the different nucleon-nucleon in- 
teraction models, this will cause a considerable theoretical uncertainty in the 
binding energies of the two-baryon systems when calculated with different 
interaction models. 

The weakening of the short range repulsion is illustrated for the two non-local 
Nijmegen potentials Nijm93 and Nijml, and for the two local potentials Nijmll 
and AV18 in Fig. 1. In the figure the matrix elements for the isospin state of 
the central interaction potential component i>i(r) — 3v TT (r) is plotted both for 
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S (CC ) - N 
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Table 1 



Quark model scaling factors for the interaction operators for two-baryon states 
formed with S c , S( cc ) hyperons and nucleons (N). The superscript + indicates lack 
of an unambiguous scaling factor. For these cases the largest scaling factor was 
chosen. 



the nucleon-nucleon and the S cc — N systems. In the latter the quark model 
rescalings in Table 1 have been taken into account. In the case of the Nijm93 
interaction model the short range attraction in the nucleon-nucleon system 
is replaced by short range repulsion in the S cc — iV system. In the case of 
the Nijml interaction model this potential matrix element becomes entirely 
attractive in the case of the S cc — N system, while in the case of the two (local) 
potentials Nijmll and AV18 the short range repulsion is almost wiped out. 

The two-baryon states of E,' c and S cc baryons differ from the deuteron in that 
the long range pion exchange interaction is also weaker by more than an order 
of magnitude. Due to this weakening this interaction is not the main source 
of binding, as it is in the case of the deuteron. 



5 




Nijml 



2 



AV18 



1 2 1 2 

r [fm] r [fm] 

Fig. 1. Central components of the potentials (v c = v i ■ 1 + v TT fi ■ fj) for isospin 
0. Solid and dashed lines show the deuteron potentials and the potentials rescaled 
with the quark model scaling factors for the H cc — N system, respectively. 



In Table 2 the binding energies for isospin combinations of E' c — N and E cc —N 
states — calculated by solving the Schrodinger equation with the rescaled 
versions of the three Nijmegen [5] and AV18 [6] models for the nucleon-nucleon 
interaction (without the isospin breaking electromagnetic components, which 
depend on the baryon charge) — are listed. For comparison, the corresponding 
binding energies that are obtained with the rescaled class of AVn' interactions 
are also listed in the table. To obtain these results only the baryon mass 
independent scaling factors from Table 1 were used. 

The very large binding energies obtained with the rescaled AV18 interaction 
models are notable. These arise from the strong L 2 interaction component in 
the AVI 8 potential. 

Due to the strongly attractive central component of the rescaled Nijml in- 
teraction, unrealistic deeply bound states with binding energies in the range 
of 1 — 10 GeV are found with that interaction model. For some of the other 
systems, that will be discussed later, the rescaled Nijm93 interaction shows a 
similar behavior. In the tables such results will be denoted by the entry * * *. 

In Table 3 the corresponding calculated binding energies for two-baryon states 
with isospin of the form S — S cc and S c — S cc as well as H cc — are listed. 
Here the mass of the was taken to be 9 GeV. The rescaled AV18 potential 
implies substantial binding energies also in the case of these states. Concerning 
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Nijmll 
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-17.2 
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— 


-1.0 


AV6' 
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— 


AV4' 


-3.0 


— 


AVX' 
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AV2' 


-0.3 


-0.6 


AVI' 







Table 2 

Binding energies of H c — A-type states with I to t = calculated with different 
interaction models and the quark model scaling factors in Table 1. The entry *** 
indicates unrealistic deeply bound states. 

the results obtained with the Nijmegen potentials, it can be seen that the two 
non-local versions, Nijm93 and Nijml, do not predict bound states and in most 
cases produce unrealistic results. The local version, Nijmll, gives rise to bound 
states. These binding energies are considerably smaller than those obtained 
with the AV18 interaction. 

The substantial spread in the calculated binding energies is a direct conse- 
quence of the large differences in the short range parts of the different po- 
tential models, which is accentuated in the case of the S— states. In Figs. 2 
and 3 the diagonal S— and D— state potentials for the S cc — N systems are 
shown for the different potential models. The Nijml potential model has a 
very attractive D— state interaction, which leads to the unrealistically large 
calculated binding energies. A comparison to the results that are obtained 
with the earlier Paris potential model [10] emphasizes this point. This poten- 
tial model leads to considerably smaller values for the binding energies: no 
binding for the the S cc — S cc and the S cc — systems and only 2 MeV for 
the — ^bb system. That the substantial cancellation between the strongly 
attractive and repulsive components in the Nijml potential model may lead 
to peculiar results in the extension to hyperons has been noted before [11]. 

The large binding energies that are associated with the L 2 component of the 
rescaled AV18 potential suggest that the quark model scaling factors in Table 1 
may be oversimplified in that they do not take into account the possibility of 
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Table 3 



Binding energies of two baryon states of S, H^. and H cc hyperons with I t ot = as 
obtained with the quark model scaling factors in Table 1. The superscript * indicates 
that more than one bound state exists. 




l 2 l 2 

r [fm] r [fm] 



Fig. 2. Full potentials in the S— wave channel (t>oo)- Solid and dashed lines corre- 
spond to the deuteron potentials and the rescaled potentials for the E cc — N system 
in the isospin state, respectively. 
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1 2 1 2 

r [fm] r [fm] 



Fig. 3. Full potentials in the D— wave channel (^22)- Solid and dashed lines depict 
the deuteron potentials and the rescaled potentials for the H cc — N system with 
isospin 0, respectively. 



baryon mass dependence in the scaling factors of the interaction components, 
which depend explicitly on angular momentum. The scaling factors in Table 1 
are obtained algebraically by application of the quark model wave functions. 
They therefore do not take any account of the fact that in dynamical models 
of the nucleon-nucleon interaction, the components which contain the angular 
momentum operator L, always contain an inverse power of the nucleon mass. 
This dependence on an inverse power of the nucleon mass plays a crucial role 
in the determination of the scaling of the nuclear interaction in the large N c 
limit [12,13]. 

In a quark model based derivation of the nucleon-nucleon interaction, it will, 
however, be the constituent quark and not the baryon masses which arise in 
the angular momentum dependent components. As in the present calculation, 
where only the interactions between the light flavor quarks are considered, 
the scaling factors in the table do not depend on the baryon mass. In order to 
investigate the numerical sensitivity to this issue, the possibility of a baryon 
mass dependence for the scaling factors, which are associated with the angular 
momentum dependent operators in (4), was investigated. 

For the L ■ S, L 2 and (L ■ S) 2 potential components for baryons of unequal 
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mass we consider the additional mass dependent scaling factor: 

2" 



1 + 



m 

m) 



(6) 



and for interactions that involve the quadratic spin-orbit interaction operator 
Q12 the mass dependent scaling factor is taken to be: 

-V • (7) 



Here m is the mass of the light baryon and M is the mass of the heavy baryon. 
The scaling factor 6 is that for a scalar interaction, and is employed here only 
for a qualitative estimate. 

In order to study the significance of this issue, the binding energies of the 
N — E' c and iV — S cc systems, as obtained when the scaling factors for those 
components are allowed to depend on mass as in (6) and (7), are shown in 
Table 4. With mass dependent scaling the AV18 potential does not lead to any 
bound state. In contrast the Nijmll interaction model result for the binding 
energies is only slightly smaller when the mass dependence of the scaling 
factors is taken into account. This is a consequence of the fact that the angular 
momentum dependence of the Nijmll interaction model is weak in comparison 
to that of the AV18 interaction. 

In the case of the isospin two-baryon states of the form S — S cc , E' c — 
S cc and S cc — S&& the results are similar, as shown in Table 5. The binding 
energies calculated with the AV18 are considerably smaller, when the mass 
dependence of the scaling factors is taken into account. This result indicates 
that more confidence should be given to the results that are obtained with the 
Nijmll potential model, which has a weaker dependence on angular momentum 
than the AV18 interaction. This conclusion is also supported by the fact that 
the angular momentum independent AVn' potentials give rise to far smaller 
binding energies than the AV18 interaction. The results for the two non-local 
Nijmegen potentials Nijm93 and Nijml remain unrealistic even when the mass 
dependent scaling is considered. This seems to indicate that these potentials 
are not very well suited for the rescaling procedure used here. 

These numerical results suggest that deuteron-like weakly bound states are 
formed by nucleons and both S c and E cc hyperons. Somewhat more strongly 
bound deuteron-like states are likely to be formed by S c and S cc hyperons 
and such bound states are also likely to be formed by S hyperons and the 
charm S c and double-charm S cc hyperons. The corresponding states formed 
with the beauty hyperons E b and (yet to be discovered) S&j, will be deeply 
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Table 4 

Binding energies of H c — N and H cc — N with I tot = as obtained with the quark 
model scalings in Table 1 with the additional mass factors in Eqs. (6) and (7) in 
the angular momentum dependent terms. 

bound. As noted above, the Nijmll interaction should be expected to give the 
most realistic estimates for the binding energies. 

Above, only the diagonal interactions in the two-baryon states formed by 
heavy hyperons have been included. Transitions between different combina- 
tions of charm hyperons are however possible through charm exchange interac- 
tions. Because of their short range they will be much weaker than the diagonal 
interactions that are mediated by light flavor exchange mechanisms. 

Two-baryon states formed by S cc may couple to A c — Q ccc states. If the total 
mass of the latter is lower than the mass of the bound S cc — S cc state, the 
latter state will only be metastable. While the mass of the Vl ccc is unknown 
experimentally, several quark models suggest that the mass of the A c — Vt ccc 
state may be lower by some 130-160 MeV [14,15]. In this case the AV18 inter- 
action, without mass dependent scalings, leads to S cc — S cc bound states but 
the Nijmll interaction only to a metastable state. 

In the case of the isospin state of the S cc — N system the charm exchange 
interactions couple it to the A+ — A+ system, which has a mass that is ~ 200 
MeV higher than that of the S cc — N bound states. Consequently the predicted 
S cc — A^ bound states are bound rather than metastable. 

The isospin state of the E' c — N system is also bound. It is coupled to the 
£ — S c system. The states of this system have ~ 130 MeV higher masses than 
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Table 5 

The calculated binding energies of E — E cc , H c — H cc and H cc — E&& states with /tot = 
a obtained with the quark model scaling factors in Table 1 with the additional mass 
factors in Eqs. (6) and (7) in the angular momentum dependent terms. 

the E' c — N states. 

In contrast the calculated "bound" states of two E! c hyperons, which may 
couple to the S — S cc state by charm exchange will be only metastable, as 
the mass of the latter is lower by ~ 370 MeV. The states of the latter are 
consequently bound and not metastable. 



3 Two-baryon states with isospin 1 charm hyperons 

3.1 Baryon-baryon interactions with isospin 1 baryons 

The scaling approach that is presented here can be readily extended to isospin 
1 baryons. As long as only baryons with spin 1/2 are considered it is only the 
isospin dependent part of the model that has to be modified. 

The quark model scaling factors can be calculated as in the previous section. 
The scaling factors for the isospin dependent interaction components refer to 
the scaling relative to the sum over the isospins of the constituent quarks, 
J2q=i T q and not relative to the baryon isospin operators r or T. Therefore, 
the procedure applies independently of the isospin of the considered baryons. 
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The isospin 1 baryon operators T 1,2,3 satisfy the SU(2) algebra [T p ,T q ] = 
ie pqr T r , without any factor \ on the right hand side as in the case of the 

Pauli isospin matrices ([|t p , \r q \ = ie vqr \r r ). Hence, the isospin operator T 

— * 

for baryons with isospin 1 corresponds to 2T = J2 q ^ q - The factor of 2 on the 
l.h.s. is included in the scaling factors employed here. This is important when 
the NN potentials are be expressed in terms of baryon operators. The isospin 
| operators Tj ■ Tj must then be replaced by Tj • ?j or Tj • Tj, respectively. The 
"missing" factors 2 or 4 are generated by the scaling factors. 

The scalings of the matrix elements for the light quark operators between £+ + 
or E° hyperons and nucleons in the quark model are: 



(E + + ,o |El , |Sc++ ,o ) = 2 (Ar|El9|iV)) 

q J q 

<£ c ++ '°|£^|£ c ++ '°> = \{N\Y,*l\N), 

q ^ q 

(S c ++ '°lE^|S c ++ '°) = 2(iV|5:r/|iV), 
q q 

< K +fi I E *W I s c ++ -° ) = U N \ E °M I n). (8) 

q q 

The isospin-z projection of the £+ hyperon is zero. Thus, the isospin dependent 
scalings are different for this charge state: 



6 q 





Eil 


|£ c + > = 




E-f 

<? 


|e+> = 




E^f 

<? 


|E+> = 






r/ 1 E c + 



(9) 



Binding energies will be calculated for isospin eigenstates that are constructed 
from the charge states. Therefore, non-diagonal matrix elements also have to 
be considered. The scalings for the matrix elements between and £+ 

hyperons and matrix elements between proton and neutron states are: 



(£r'°iE^ c + >=±^<piE^>, 



(s c ++ '°iE^/ls c + ) 



2V2 



p I E a i T ! I n > 



(10) 
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/,From the differences in the isospin dependent scalings in Eqs. (8)-(10) it 
follows that the full scaling behavior of isospin eigenstates cannot be described 
by the quark model scaling factors alone. The crucial point is that the structure 
of the isospin eigenstates changes when isospin 1/2 baryons are replaced by 
isospin 1 baryons. There are then three E+, S°) charge states instead 

of two (p, n). This change of the structure — which manifests in different 
Clebsch-Gordan coefficients — must be seen as a part of the scaling procedure. 
The quark model scaling factors are derived for the individual charge states. 
They have no information about the arrangement of the charge states into 
isospin eigenstates. The total scaling will thus be a combination of the quark 
model scaling factors and the Clebsch-Gordan coefficients, which has to be 
determined individually for each isospin eigenstate. 

In the following sections it is shown that the full scaling of a isospin depen- 
dent interaction component can be split up into a quark model scaling factor 
and another scaling factor that is related to the explicit value of the matrix 
element of the two-baryon state considered. The matrix element relations for 
the relevant isospin eigenstates are: 



( 1; ±m | Ti ■ Tj | 1; ±m) = 1 , m — 0, 1 , 

( 0; | n ■ Tj I 0; ) = -3 , (11) 
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(-;m\T i -T j \-;m} = l, 



m — -, - 



2' 1 1 Jl 2' ' ' 2'2' 

i ± _ \ -L i ' Tj -: ±- 

2 2 1 3 1 2 2 



1 ;±i|r i -7*-|i;±i> = -2, (12) 



{2\±m\Ti'Tj\2\±m) = l, m = 0,1,2, 
( 1; ±m | Tj ■ Tj | 1; ±m ) = — 1 , m = 0, 1 , 

(0;0|f l -f i |0;0) = -2. (13) 



Note that the composition of the eigenstates could be ignored when only 
isospin 1/2 baryons were considered. The isospin eigenstates of nucleons and 
2(cc,fefe) baryons are constructed identically, i.e., using the same Clebsch-Gordan 
coefficients. Therefore, the full scaling information is contained in the quark 
model scaling factors. 
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3.2 Bound states of two S c hyperons 



The E+ + and the S° hyperons form isospin symmetric two-body states with 
isospin 2 and one antisymmetric state with isospin 1. These are explicitly: 



2;2) = 
1,0) = 



1 

71 



>, I 2, — 2 > = 



(14) 



The states | 2; ), | 2; ±1 ), | 0; ) and 1 1; ±1 ) are linear combinations of all 
the S c charge states. Their explicit expressions are: 



2;0) = 

2;i> = 
2;-i) 



l 
1 

71 



1 

71 



s c ++ s°) + 2|s c + s+) + |s°s c ++ ; 

) + I S c + S+ + ) 



(15) 
(16) 
(17) 



and 



lil>= ^[ 






1; - 1>= ^ 









) + I e^e c 



(18) 
(19) 
(20) 



respectively The isospin independent scaling factors for those states, which 
can be derived directly from Eqs. (8), are listed in Table 6. The full scaling 
behavior of the isospin dependent interaction components is found by calcu- 
lating the scaling between each of the two-baryon state matrix elements and 
the corresponding isospin symmetric or antisymmetric two-nucleon matrix el- 
ement. 

The simplest case is the | 2; 2 ) state, since the structure of this isospin eigen- 
state is identical to the structure of the | 1; 1) state that is formed of two 
isospin 1/2 baryons. Using Eqs. (8) it is found: 



(2;2|E^-E^|2;2) = (s c ++ |E^I^ + )-(s 6 ++ |^ 

vi q q 

= 2<p|£**|p>-2<p|E**IP> 



r<ME+ + 



15 



=4<i;i|E*f-E*?|i;i> 



NN 



[21) 



The state | 0; ) mixes the charge states. Thus, all the scalings from Eqs. (8)- 
(10) enter into the two-body matrix elements: 



(0;0|£tT-£t7|o ; o; 



q q 
q q 

-(s c + iE^I^ + )-(s c t|E^|s°) 

q q 
2(p|E^IP>-2< n lE^l n > 



-2V2(p\ J2^\n) ■V2{n\J2^ Q \p] 

q q 

=4-|(o ; oiE^f-E^lo;oW- ( 22 ) 

<3 p g 



This calculation can be repeated analogously for all the states (14)-(20). It is 
found that the scaling only depends on the total isospin of each state and not 
on the isospin-z projection. The explicit results for the scaling factors are: 



(2;±m|E^-E^|2;±m)=4.1(l;l|E^-E^I 1 ; 1 ) 



NN ■ 



(l;±m| E^T-E 7 ? |l;±m>=4- -<0;0| E*f • E*i I °5 °>™ ■ 

p q p q 



(0;0|E^-E^1lO;0) = 4.-(0;0|E^f-E^|0;o: 



NN ■ 



(23) 



The factor 4 can be interpreted as the quark model scaling factor that is 
derived from the matrix element scaling in (8). The meaning of the factors 1, 
| and | becomes clear from Eqs. (11)-(13). These isospin dependent factors 
are given by the ratio of the matrix elements for the considered two-baryon 
state and the corresponding symmetric or antisymmetric two-nucleon state: 



( 2; ±m | Tj • T 3 


2;±m) 


(i;i| 


n ■ ^ 


l; 1 /nn 



m = 0,1,2, 



(24) 



( 1; ±m 


1 T ■ T 


1; ±m ) 


(0;0| 


Ti ■ Tj 


0;0) NN 



m = 0, 1 



'3 



(0;0 


1 T ■ T 


0;0) 


(0;0| 


Ti ■ Tj 


0;0) NN 
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S c -iV 


— 


Operator 


Scaling factor 


1 


2/3 


4/9 


n • fj 


2 


4 


G{ - Gj 


4/3 


16/9 


{Oi-Oj){fi-fj) 


4/5 


16/25 


Sij 


4/3 


16/9 


Sijif • fj) 


4/5 


16/25 


L-S 


4/3+ 


8/9 


L-S(fi-fj) 


2+ 


8/5 


I? 


2/3 


4/9 


L 2 (n-f,) 


2 


4 


L 2 (ai ■ Gj) 


4/3 


16/9 


L 2 (ai ■ Gj)(fi ■ fj) 


4/5 


16/25 


(L • S) 2 


4/3+ 


16/9+ 


{L.Sf{f.fj) 


2+ 


4+ 


Ql2 


4/3 


16/9 


Qi2(n ■ fj) 


4/5 


16/25 



Table 6 

Quark model scaling factors for the interaction operators for the two-baryon states 
T, c — N and S c — S c . The isospin dependent scaling factors must be multiplied with 
the appropriate factor from Eqs. (24) or (31) to obtain the full scaling factors. 

For the other isospin dependent interaction components, like {B i -Bj){f i -fj), the 
same behavior is found by a similar calculation. The complete scaling factors 
are products of the scaling factors that can be derived from Eqs. (8) and the 
isospin dependent scaling factors (24). The quark model scaling factors are 
listed in Table 6. The scaling factors that could be calculated from Eqs. (9) 
and (10) do not appear explicitly in the total scalings. They have merged with 
the Clebsch-Gordan coefficients to provide an unabiguous scaling behavior of 
the isospin eigenstates. 

Note that the S+ + — S+ + and S° — S° light flavor quark interactions scale 
from the corresponding nucleon-nucleon interactions in the same way. The 
flavor breaking terms are insignificant as S+ + and £° have almost the same 
mass (Am = 0.4 MeV). Hence, the results for bound states with different 
isospin-z projections will be the same. 
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|2,m) 


1 1, m) 


1 o,o> 


Potential 


Binding 


Energy [MeV] 


Nijm93 


-66.6 




#*# 


Nijml 


— 


*** 


*#* 


Nijmll 


— 


-53.7 


-285.8 


AV18 


-41.1 




#*# 


AV8' 


— 




-16.1 


AV6' 


— 




-10.8 


AV4' 


— 


-7.3 


-87.4 


AVX' 


— 


-2.8 


-53.3 


AV2' 




-8.3 


-58.5 


AVI' 


-0.7 


-0.7 


-0.7 



Table 7 

Calculated binding energy for S c — S c states obtained with the quark model scalings 
of the interaction operators. 



The calculated binding energies for the two-charm hyperon states are listed 
in Table 7. For the symmetric isospin 2 states the Nijm93 interaction and the 
AV18 interaction lead to bound states, while the Nijmll and most of the AVn' 
potential models yield no bound states. In contrast, the Nijmll and many of 
the AVn' potentials give rise to bound states for the antisymmetric isospin 1 
and states. The isospin state is much stronger bound than the isospin 1 
states. While the binding energies are rather high there is no convergence on 
the value of the binding energy for the different potentials. 

The bound states of two S c hyperons found in Table 7 for some of the inter- 
actions are only metastable. The isospin 2 states can couple to the lower lying 
S cc — A(1232) states that have the same isospin by the short range charm 
exchange interaction. In the same way the isospin 1 and states can couple 
to the even lower lying S cc — N states. Even the the strong binding for the 
isospin states from the AV18 interaction cannot compensate the large mass 
difference of ~ 500 MeV. 

In contrast the bound states of S c hyperons and nucleons that will be discussed 
in the following section are stable. They can only couple to the A c — A(1232) 
and S c — A(1232) states which have ~ 120-300 MeV higher masses. 
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3.3 S c — N-type bound states 



The simplest combinations of nucleons and E c hyperons are the two symmetric 
isospin 3/2 combinations: 



3 3 
2' 2 



3 
2' 



(25) 



The isospin 3/2 states with 1/2 and —1/2 isospin-z projections are the two 
symmetric linear combinations: 



-■-) 


1 r 


2' 2 / 


~7^ 


3 1 




2'~2 


) = 7z 



S++n) + v / 2|S+p) 



y/2\ E+n) + |E°p> 



In addition there are the two antisymmetric isospin 1/2 combinations: 



(26) 
(27) 



1 1 

2' 2 
1 1 
2' ~2 



1 

7! 



V2|S 



++ 



n , 



1 

71 



E+n)-x/2|E°p) 



(28) 
(29) 



As in the previous section, the isospin independent scaling factors can be de- 
rived directly from Eqs. (8). The full scaling behavior of the isospin dependent 
interaction components has to be determined by calculating the scaling be- 
tween each of the two-body matrix elements and the corresponding isospin 
symmetric or antisymmetric two-nucleon matrix element. Again, the scalings 
will not depend on the isospin-z projections but only on the total isospin of 
each state. From calculations similar to those in Eqs. (21) and (22) full scaling 
factors are found: 



(|;±mlE^-E^l|; ± "»> = 2 - 1 < 1 ; 1 IE^-E^|i;iW. 

p <i p q 

(^±^IE^-E^I^±^) = 2-|(o;oiE^-E^lo;o)^.(3o) 

Z p q ^ " p q 

The factor 2 on the r.h.s. of these equations is the quark model scaling factor 
derived from the matrix element scaling in Eqs. (8). Analogously to Eqs. (24) 
the isospin dependent factors 1 and | are the ratios of the matrix elements 
for the considered two-baryon states and the corresponding symmetric or an- 
tisymmetric two-nucleon state: 
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|3/2,m 


|l/2,m) 


Potential 


Binding 


; Energy [MeV] 


Nijm93 


-5.3 


— 


Nijml 


— 


-227.3 


Nijmll 


— 


-16.9 


AV18 


— 


— 


AV8' 


— 


— 


AV6' 


— 


— 


AV4' 


— 


-4.9 


AVX' 


— 


-1.9 


AV2' 




-4.3 


AVI' 


-0.2 


-0.2 



Table 8 

Calculated binding energies for T, c —N states as obtained with quark model scalings. 



1 = 



( §; ±m 


1 T i ■ *j 


|; ±m) 


< i; 1 1 

\ 2' 2 


n ■ Tj 
\fi-f 


l 

i 


; 1 )nn 

2' 2 / 


(0;0| 


n ■ 7j 


0;0)^ ' 



(31) 



The scaling factors for the other isospin dependent components of the inter- 
action are found in the same way. The quark model scaling factors for all 
interaction components are listed in Table 6. As discussed for the two-charm 
hyperon states the results for the binding energies will be identical for states 
with different isospin-z projections. The same scaling factors are used and the 
small mass differences between the £+ + and TP C are ignored here. 

In Table 8 the calculated binding energies are listed. Only the Nijm93 and the 
AVI' potentials suggest the existence of very weakly bound isospin 3/2 states. 
Two of the Nijmegen potentials and some of the AVn' potentials suggest the 
existence of bound isospin 1/2 states. The Nijml potential predicts a much 
higher binding energy than the other interactions. For the Nijm93 and the 
AV18 potentials no bound states are found. It should be noted that — as for 
the S c — E c systems — the Nijmll and the AVn' potentials prefer the existence 
of bound states with low total isospin. 
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4 Discussion 



The investigation above considered the use of well known phenomenological 
nucleon-nucleon interaction models to approximately describe the interaction 
between nucleons and heavy flavor hyperons and to explore the possible ex- 
istence of deuteron-like bound states of nucleons and heavy flavor hyperons. 
The method relies on quark model scaling factors for scaling the strengths of 
the different interaction components to the appropriate number of light flavor 
quarks in the heavy flavor hyperons considered. While straightforward in exe- 
cution the method proved to be only of qualitative value because the large vari- 
ation in the short range components of the different modern nucleon-nucleon 
interactions led to considerable scatter in the calculated binding energies of 
the two-baryon systems considered. 

Another limitation is the fact that the quark model scaling factors for the 
angular momentum dependent interaction components, which relate to the 
light flavor quarks, do not reveal the expected dependence on the hyperon 
mass, which is suggested by the large N c limit of QCD [12,13]. Because of 
this the binding energy estimates that are obtained with the Nijmegen and 
the AVn interaction models should be viewed as more reliable than those 
obtained with the AV18 interaction, which has a strong quadratic angular 
momentum interaction. 



With these provisos, the present results suggest that nucleons form bound 
states with S c hyperons in the isospin state that have binding energies in 
the range between 3 and 14 MeV. With somewhat less confidence nucleons 
are expected to form bound states with E cc hyperons, with binding energies 
between 1 and 34 MeV. Such deuteron-like bound states of S cc hyperons, and 
of E bb as well as of S cc and S b6 hyperons are also very likely, with binding 
energies in the range 90 - 120 MeV. Bound states of nucleons and S c hyperons 
are also likely, although their binding energy cannot be estimated with much 
confidence by the present method. Most likely states with low total isospin 
will be most strongly bound. 

Here the interactions between strange and light flavor quarks have been ne- 
glected. A more realistic description of the states that contain strange quarks, 
e.g. the S or S c hyperons, should be possible by using a rescaled baryon- 
baryon potential like the Nijmegen NSC97 potential [16]. However, due to the 
limited experimental data that is available to fit such potentials, this would 
also introduce additional uncertainties. 
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A Quark model wave functions 



The quark model wave functions of the baryons considered in this work can 
be written as: 



= V s (r)[<Kflavor)x(spin)U a (color) , (A.l) 

where the subscript s denotes a symmetric and the subscript a an anti- 
symmetric component of the wave function. 

Nucleons and S( cc ), E^, S c hyperons have spin and flavor wave functions with 
mixed symmetry. The symmetric spin-flavor part of the total wave function is 
constructed as follows: 



[4>X]s = [(pmsXms + (pmaXma] • (A.2) 

The explicit forms of spin and flavor wave functions may be constructed by 
the methods in ref. [17]: 



T ) ms = TTI + TIT -2 ITT ) , I T )r 

)n» = -^|UT + m-2TU), ll)r 



-j=\ TTI - TIT 
^2 1 ITI IIT 



(A.3) 



\p)ms = \ uud + udu — 2duu) , \p) 



ddu + dud — 2udd ) , \n\ 



—= | uud — udu ) , 
v2 



V2 



dud — ddu ) , 



(A.4) 



S+ + ) ms = - | ecu + cue - 2ucc ) , I S+ + ) ma = I cue — ecu ) 
I 5+ ) ms = --^=| ccd + ede- 2dcc) , | S+ ) ma = -J=| ede - ccd) , 



(A.5) 
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I S c + ) ma = ^_ | {ucs + scu) + {use + sue) — 2(csu + cus) ) , 

^'c + )mo = ^| ( M SC + SWc) - (SCW + WCs) ) , 

I S c ° ) ms = — | (scd + des) + (sdc + dsc) — 2(cds + csd) ) , 
vl2 

„'0 X 1 

2 



c /ma = — | (sdc + dsc) — (des + scd) ) , (A. 6) 



I 2+ )ms — 7 — I {ucd + dcu) + (udc + d-uc) — 2(cdw + cud) ) , 
V 12 

I E+ ) ma = -| (udc + due) — (dcu + -ucd) ) , 
I A+ ) ms = -| {udc — due) + {ucd — dcu) ) , 

I A+ ) ma = — | {ucd — dcu) + {due — udc) — 2{cdu — cud) ) , (A. 7) 

vl2 



1 1 

I )ms — —/= \ ddc + ded — 2cdd ) , | E c ) ma = —j= \ ddc — ded ) , 
v6 v2 

I ) ms = —= | -uwc + ucu — 2cuw ) , | £+ + ) ma = — = | uuc — ucu ) 
y6 v2 



(A.8) 
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